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& ' Abstract 



The analytic signal method via the Hilbert transform is a key tool in signal analysis and process- 
ing, especially in the time-frquency analysis. Imaging and other applications to multidimensional 
signals call for extension of the method to higher dimensions. We justify the usage of partial 
Hilbert transforms to define multidimensional analytic signals from both engineering and mathe- 
matical perspectives. The important associated Bedrosian identity T(fg) = fTg for partial Hilbert 
transforms T are then studied. Characterizations and several necessity theorems are established. 
We also make use of the identity to construct basis functions for the time-frequency analysis. 
Keywords: multidimensional analytic signals, partial Hilbert transforms, the Bedrosian identity. 
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^ ■ 1 Introduction 

Signals are carrier of information. In many applications, feature information of a signal usually persists 
£NJ ■ in the time- frequency domain [9] . The analytic signal method [10] is a classical way of defining without 

ambiguity the local amplitude and frequency of one-dimensional signals. Among others, it has been 
proven useful in meteorological and atmospheric applications, ocean engineering, structural science, 
and imaging processing, [13j [151 CES1 US EI] ■ Especially, it motivates the widely used empirical mode 
$_i . decomposition and the Hilbert-Huang transform [131 111 E] • 

The analytic signal method makes use of the Hilbert transform H defined for functions / G L 2 (R) 

as 

(Hf)(x) :=p.v~ / I^-dy:= lim - I -^Ldy, x£l, (1.1) 



KjRX-y E ->0+ 7T J\y- X \> £ x - y 

Let H 1 (R) be the Banach space of all the functions 4> G L 1 (M) such that H(j) G L 1 (M) equipped with 
the norm 

IHIhi(»o := IHIlmk) + II^IIl^r)- 

The Hilbert transform Hf of / G L°°(R) is defined as a BMO function that is in the dual of H 1 (R) by 



/ {Hf)4>dx = - [ fHcpdx, 
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The analytic signal Af of a real / G L 2 (R) is formed by adding to it its Hilbert transform as the 
imaginary part: 

Af:=f + iHf. (1.2) 

With the amplitude-phase decomposition 

{Af)(t)=p(t)e»®, t€R, 

the p{t) and 0'(t) respectively are taken as the instantaneous amplitude and frequency of the signal / 
at time t. 

The Hilbert transform in (|1.2|) has many favorable mathematical properties that account for the 
usefulness of the analytic signal method in the time-frequency analysis. Here we mention four of them. 
Let N be the set of positive integers and Z + := N U {0}. For a fixed d G N, we denote by || • || and 
(•, •) the standard Euclidean norm and inner product on R d . The Schwartz class S(R d ) consists of 
infinitely differentiable functions ip on M. d such that for all p G Z + and q G 

sup{(l + \\x\\ 2 ) p \f^(x)\ :x£R d }< +00. 

The Fourier transform (p of <p G is defined as 

■= [ ip(x)e- l ^dx, £ G R d , 

which is again a function in 5(]R a! ). The Fourier transform can be extended to the space S'(M. d ) of 
temperate distributions on R d by a duality principle [llj . 

The Hilbert transform has an equivalent definition via the Fourier multiplier — i sgn , where sgn (£) 
takes value —1,0, 1 for ^ < 0, ^ = and ^ > 0, respectively. Specifically, we have for all / G L 2 (IR) 

(Hf)X0 = -isgn(0f(0, (1-3) 
There are two important consequences of this. First of all, we have for all / G L 2 (M) that 

supp(/ + ifr/)CR+ (1.4) 

where := [0, +00). In other words, the analytic signal suppresses all the negative frequency 
components of the original signal. This justifies the physical soundness of the analytic signal method. 
Secondly, the multiplier definition (11.31) makes it possible to develop fast algorithms making use of the 
lighting fast FFT to compute the analytic signal [18]. In addition to the above two properties, it was 
shown in |28|. [29] that the Hilbert transform is the only continuous and homogeneous operator L such 
that 

L(cos(ojot + <^>o)) = sin(cjoi + 4>o), for all > 0, ^ 6 M. 

This implies that the analytic signal method is the only one that satisfies some reasonable physical 
requirements (28, 29]. Finally, the Hilbert transform satisfies the Bedrosian theorem [T]: If /, g G L 2 (M) 
satisfy either supp/ C R + , supp^ C R + or supp/ C [—a, a], supp^ C (— 00,— a] U [a, 00) for some 
positive number a, then there holds the Bedrosian identity 

[H(fg)](x)=f(x)(Hg)(x), x G R. (1.5) 

The Bedrosian identity simplifies the calculation of the Hilbert transform of a product of functions, 
helps understand the instantaneous amplitude and frequency of signals, and provides a method of 
constructing basic signals in the time- frequency analysis [H El HI EJ HHJ EH [231 E3 [36] . It has attracted 
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much interest from the mathematical community O [8l [22l [231 EOl ET], [33l [3¥l [35l [36] . Here, we mention 
an observation in [30]. It states that the Hilbert transform is essentially the only operator that satisfies 
the Bedrosian identity. 

We conclude that the analytic signal (|1 .2j) is justified by the above mathematical properties. This 
paper is motivated by the need [T2"] 152] of defining multidimensional analytic signals for the time- 
frequency analysis of multidimensional signals. Naturally, we are inclined to define the analytic signal 
of / G L 2 (R d ) through a fixed operator T : L 2 (R d ) -> L 2 (R d ) as / + Tf. Thus the definition reduces 
to the choice to the operator T. In references [5JQ2]; compositions of the partial Hilbert transforms 
are used. One of our purposes is to show that linear combinations of compositions of the partial 
Hilbert transforms are indeed the only choice from the viewpoint of similar mathematical properties 
introduced in the last paragraph. The justifications will be carried out in the next section. In Section 
3, we investigate the multidimensional Bedrosian identity T(fg) = fTg for f,g€ L 2 (R d ), where T is 
a linear combination of compositions of the partial Hilbert transforms. In particular, a necessary and 
sufficient condition and the Bedrosian theorem will be established. The necessity of the Bedrosian 
theorem will be discussed as well. In Section 4, we construct basis multidimensional analytic signals 
by the results on the Bedrosian identity. 



2 Multidimensional Analytic Signals 

Set 7L n := {0, 1, . . . , n — 1} and N ra := {1,2,..., n}, n G N. The partial Hilbert transform Hj, j G 
is defined for / G L 2 (R d ) as 



(Hif)(x) :=p.v.- [ J™ dy:= lim - f JM-dy, x G F d 

KjRdXj-yj £ ->0+ 7T J\ yj - X:j \>e Xj ~ Vj 



It is known that Hj is a bounded linear operator from L 2 (R d ) to L 2 (R d ), |26] , Moreover, there holds 

(F i /r(0 = -^gn(^)/(0, ^i d , je% (2.1) 

The purpose of the this section is to justify the using of linear combinations of compositions of the 
partial Hilbert transforms to define the multidimensional analytic signal. For simplicity, we denote 
by Hq the identity operator on L 2 (R d ). A linear combination T of compositions of the partial Hilbert 
transforms has the form 

where H a := TljeNd an d c a are complex constants. We shall argue from four points of view that 
the multidimensional analytic signal Af of / G L 2 (R d ) should be defined as 

Af:=f + Tf (2.3) 

via an operator T of the form (12. 2\\ . 



2.1 Concentration of the Frequency 

By (|1.4p . the negative Fourier frequency of a one-dimensional analytic signal is suppressed while the 
energy of the positive frequency is doubled. For some applications, it is desirable to restrict the 
frequency to a certain d-hyperoctant, especial the first one, [12J. Following the notations of [30j, we 
let u k , k G N 2 d, be the extreme points of the cube [—1, l] d . Assume that u 1 and u 2 are the point each of 
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whose component is 1 and — 1, respectively. The Euclidean space R is divided into 2 d-hyperoctants 
as 

R d = \J Q k , 

k€N 2d 

where Qk := {£ : £ G R rf , > 0, j G N^}. Suppose that we have chosen the fc-th <i-hyperoctant and 
desire that for each / G L 2 (R d ) its multidimensional analytic signal Af be such that 

wm) = { 2 '{«>- o « h ^ (2.4) 

The above equation has the following equivalent form 

(A/m) = ho n ^ + "* s § n &))> £ G 

By ()2.ip , the above equations holds true for all / G L 2 (R rf ) if and only if A is defined by (|2.3p with T 
there being given as 

T= HiHo + i^H^-Ho, 

which satisfies (12.21). 



2.2 Computational Advantages 

As a consequence of (|2.ip . it can be seen that a bounded linear operator T has the form f)2.2f) if and 
only if there exists some function m G L 00 ^^) that is constant on each d-hyperoctants such that 

(Tfy=mf, f G L 2 (R d ). (2.5) 

In other words, operators (|2.2p are given by a very simple Fourier multiplier. Therefore, efficient 
numerical algorithms based on fast multidimensional discrete Fourier transforms (see [2] and the 
references cited therein) can be developed for the computation of the multidimensional analytic signal 
(|23j) . 

2.3 Physical Requirements 

Following [28, 29], we ask what continuous linear operators T preserve the class of complex sinusoids. 
To make this precise, we need to extend the definition of partial Hilbert transforms. For each a G 
we define U a : L 1 (R d ) -»• S'(R d ) by setting for each / G L x (R d ) as 

:= IT (-ir Ml ' ai) [ fHatpdx, f G 5(R d ), (2.6) 

where for (j) in a locally convex space V and «£ V* := fct(<?!>). Since % a is bounded from 5( 

to L°°(R d ), (USD is well-defined. Denote by H 1 (R d ) the Banach space of functions / G L x (R d ) such 
that Hjf G L 1 (R d ) endowed with the norm 
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The dual of M 1 (IR a! ) is the space BMO(M rf ) of locally integrable functions on R d with bounded mean 
oscillation [26], which embeds continuously into S'(R d ). Note that C BMO(lR d ). It is clear 

that n a is bounded from H 1 (R d ) to H 1 (R d ) for each a G Z^ +1 . Thus, we define for / G BMO(R d ) 
"Hq,/ again as a BMO function by 

(V,J) == II (-l) min(1 ' ai) (/>^) 5 5 G M 1 ^). (2.7) 

For £ G M d , we set sgn (£) := (sgn (£j) : j G and shall characterize operators that preserve the 
class of complex sinusoids. 

Proposition 2.1 Suppose that T is a linear operator from L co (R d ) U L 2 (R d ) to S'(R d ) that is con- 
tinuous in the topology of S'(R d ) and satisfies for some function A : C d — > C that 

Te iM = A( sgn (w))e i(av) , for all £ G R d . (2.8) 

Then T has the form 112. ty) . 

Proof: Let / G L 2 (R d ). Since § := span{e 1 ^') : uj G R d } is dense in S'(R d ), there exists a sequence 
f n G § that converges to / in iS^l^). By the continuity of T, Tf n converges to Tf. Thus, (Tf n )" —> 
(Tf)". Denote for uj G R d by 5 W the delta distribution at uj. We observe that 

(e^T = (2n) d 6 u . 

Consequently, 

(Te^-)r=(27r) d A(sgn( W ))^. (2.9) 

Let cp be an arbitrary function in 5fc(M a! ) := {p G 5(M a! ) : suppt/? C for some k G N 2 d. Then we 
see that A( sgn (•))<£ G S(R d ). Assume that /„ = EjeN fc Cnje^^. We get by that 

{(TfnY, i P)= ^2 c nj(2n) d \( sgn (uj nj )p(uj nj ) = (f n , A( sgn (•))(/?) . 

Since 

lim <(T/ n )» = ((T/)» and lim (f n , A(sgn (-)V) = (/, A( sgn (•))</?>, 
we have for all cp G 5fc(M d ), G N 2 d that 

((T/)» = (A(sgn(-))/^). 



Since the linear span of Uk^n 2d Sk(R ) is dense in S(R ), the above equation remains true for all 

cp G 5(R d ). Therefore, (Tf y = A( sgn (•))/. We complete the proof by pointing out that a linear 
operator T : L 2 (R d ) — > L 2 (R d ) has the form (|2.2p if and only if it is defined by a Fourier multiplier 
that is constant on each Qk, k G N 2 d. □ 



2.4 The Bedrosian Identity 

Considering the importance of the Bedrosian identity in the time- frequency of one-dimensional signals, 
it is natural for us to expect a linear operator T in the definition of the multidimensional analytic 
signal (|2.3p to satisfy a multidimensional Bedrosian identity. Assume that a bounded linear operator 
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T : L 2 (M. d ) — > L 2 (M. d ) is engaged to define multidimensional analytic signal by (|2.3p . First of all, we 
shall require that T is translation invariant, namely, for all x G M. d and / G L 2 (JH d ) 

T(/(--x)) = (r/)(.-x). 

The reason is that there is usually a delay in time of an input in signal processing. Thus we would 
like the processor T to have the property that a delay in the input will simply cause the same delay 
in the output. 

Now we are seeking extensions of the Bedrosian theorem that a bounded linear translation invariant 
operator could possibly satisfy. We say that a bounded linear operator T : L 2 (M. d ) — > L 2 (M. d ) satisfies 
the type one Bedrosian theorem if the Bedrosian identity 

T(fg) = fTg (2.10) 

holds true for all /, g G that satisfy 

supp/ C Y\ [— a j, a j], suppg C Y\ R\ ( — f° r some a = (dj > : j G N^). (2-11) 

jen d jen d 

The following results was proved in |30j. 

Lemma 2.2 A bounded linear translation invariant operator T : L 2 (M. d ) — > L 2 (R d ) satisfies the type 
one Bedrosian theorem if and only if it is of the form 

Thus, if we desire that the operator T satisfy the type one Bedrosian theorem then it must be a 
linear combination of the compositions of the partial Hilbert transforms. A more natural extension of 
the Bedrosian theorem is the type two Bedrosian theorem, that is, f|2. lOj) holds true for all f,g*E 5(M ) 
satisfying 

supp / C fl(0, R), supp gCR d \ B°(0, R), for some R > 0, (2.12) 

where B(xo,r) := {x G M. d : \\x — xq\\ < r} for xo G M. d and r > 0, and B°(xo,r) denotes its interior. 
When d = 1, the type one and type two Bedrosian theorems are identically the same. However, we 
shall prove that for d > 2 there exist no nontrivial operators that satisfy the type two Bedrosian 
theorem. 

Theorem 2.3 Let d > 2. Then there does not exist a nontrivial bounded linear translation invariant 
operator T : L 2 (M. d ) — > L 2 (M. d ) that satisfies the type two Bedrosian theorem. 

Proof: Let T : L 2 (M. d ) — > L 2 (IR c( ) be bounded linear and translation invariant. Assume that it satisfies 
the type two Bedrosian identity. Let /, g be arbitrary functions in 5(M rf ) that satisfies (I2.1ip . Let 
R := ||a||. Then we observe that (I2.12p holds true. By the assumption, the Bedrosian identity (|2.10l) 
holds. We conclude that (I2.10p holds for all f,g G S(R d ) satisfying (|2. 11 j) for some a := (oj > : j G 
Nd). By Lemma l2.2| T must be of the form (12. 2p . Thus, T is defined by (12.51) via a Fourier multiplier 
m G L°°(R rf ) that is constant in each d-hyperoctant. Let be the constant value that m takes in 
Q k , k G N 2 d. 

Assume that T is neither the zero operator nor a multiple of the identity operator. In other words, 
"ifci 7^ "ifc 2 f° r some k%,k2 G N 2 d. Since d > 2, there must also exist some k% G N 2 d such that 
m k\ 7^ m k 3 or m k2 7^ m k s - There hence exists a pair h,l2 G N 2 d such that ^ mi 2 and v^v^ = 1 
for some jo G N^. For notational simplicity, assume that l\ = 1 and jo = 1. Since Z 2 7^ 1, the set 
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{j G : v l ? = — 1} is nonempty. As a consequence, there exist positive constants e, r\,t2 so that 
C 1 ^ 2 G M d defined by £ x := (1 - § , n, n, . . . , r x ) and 

f 1-e, fc = l, 
I -r 2 , ^ 2 = -l- 

satisfy H^H < 1, ||£ 2 || > 1 and Hf 1 -£ 2 || < 1. Choose r > so small that B(^ -£ 2 ,r) C Qi n 5(0,1) 
and 5(£ 2 ,r) C n (M d \ 5(0, 1)). One can construct nonnegative /, <? G S(M. d ) such that supp/ C 
B^ 1 - £ 2 ,r), supp 5 C 5(£ 2 ,r), /(f 1 - £ 2 ) > 0, and g^ 2 ) > 0. Then we have that 

supp/ C B(0, 1), supp g <ZR d \ B°(0, 1). 

Since T satisfies the type two Bedrosian theorem, (12. 1QH holds. Applying the Fourier transform to 
both sides of (pTTO]) yields that 

/(£ - V)aiv)imi0 ~ min))dr) = 0, for all £ G R d , 

which has the following form at £ = 



(mi - m; 2 ) / /(£ 1 - i])g(ri)dr) = 0. 

However, since mi 7^ m; 2 , /(£ — £ 2 )<?(£ 2 ) > 0, and /(£ x — ■)# is nonnegative, 

(mi - m/ 2 ) / fiC 1 - r])gir])drj / 0. 

We hence come to a contradiction and complete the proof. □ 

Based on the above mathematical considerations, we conclude that the only operators that should 
be used in (]2.3p to define the multidimensional analytic signal should be linear combinations of the 
compositions of the partial Hilbert transforms. 



3 Multidimensional Bedrosian Identities 

The purpose of the this section is to study the Bedrosian identity fl2. lOf) for the time- frequency analysis 
of multidimensional analytic signals. By the discussion in the last section, we shall assume from now 
on that T is of the form (|2.2p . In other words, there exist constants m^ G C, k G N 2 d such that T is 
defined by (|2.5p through the Fourier multiplier m G L°°{R d ) given as 

m(0 := m k , £ G Qk, k G N 2d . (3.1) 

Denote by T* the adjoint operator of T on L 2 (K. d ), that is, 

(Tf,g) L2m = if,T*g) L * m , f,g G L 2 (R d ), 

where (•, -)i 2 (R d ) is the inner product on L 2 (M. d ). It is clear that T* is defined by the Fourier multiplier 
m. By definition f|2.6[) . as a linear combination of compositions of partial Hilbert transform, T is 
well-defined on L 1 (M d ) with range in S'(M. d ) by 

(Tf,<p) = (f,T^), V eSiR d ). (3.2) 
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We shall start the characterization of the Bedrosian identity with a lemma on a property of the 
operator T. To this end, we recall that the Fourier transform S of a temperate distribution S G 
is again a temperate distribution satisfying 

(S,<p) = {S,0), ^eS(R d ). (3.3) 

Lemma 3.1 For all f G L 1 (R d ), (Tf)"=mf. 

Proof: Let / G L 1 (M d ). We proceed by equation ([331) and <E2]) that for all ip G S{ 



(Tf,p) = (Tf,0) = (f,T*0). (3.4) 



Set /i := (T*tp)~, the inverse Fourier transform of T*ip. Note that h G -L 1 (R rf ). Thus, by ([3? 

<(T/)» = / /k= / fhdx.. 

JR d JR d 

We compute that 



h = (T*0)' = -j—- 2 (T*(py = (T* (<p)J = m<p = nvp. 
Combining the above two equations proves the lemma. □ 

Theorem 3.2 Let T be defined by the Fourier multiplier \3. Then f,g G L 2 (R d ) satisfies the 
Bedrosian identity )12.10\) if and only if 

J2 (mk-mj) f /(£ - rj^Wv = 0, £ G Qk, k G N 2 d. (3.5) 

Proof: Since the Fourier transform is injective from 5'(M d ) to 5'(M d ), (12.10P holds if and only if 

(T(fg)T=(fTgy. 

By Lemma 13.11 and writing the Fourier transform of a product of two functions as their convolution, 
the above equations is equivalent to 



m(0 / /(£ " V)9(v)dv = /(£ - Ti)m(ri)g(Ti)dn, a.e. £ G K d . (3.6) 

jR d JM d 

Since the left hand side above is continuous about £ G M. d , (|3.6[) holds if and only if 

m(0 / /(£ - = / /(£ " ri)m(ri)g(r))dri, for all £ G M d . (3.7) 

By dSH]), (J321) has the form (pT5) . The proof is complete. □ 

In the one-dimensional case, various characterizations of the Bedrosian identity have been proposed 
in the literature. For instance, the above theorem for d = 1 was proved in [36]. A similar result was 
presented in [3] under the assumptions that /, g G L 2 (M) n L°°(R). Other characterizations can be 
found in [331 ES]- 

We now turn to sufficient conditions for the Bedrosian identity. We first prove the following 
generalization of the Bedrosian theorem, which was first proved in [25J for the case when T is the total 
Hilbert transform YljeUd^' 
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Proposition 3.3 If f,g G L 2 (R d ) satisfy for some a := (a,j > : j G N^) and b := (bj > : j G N^) 

supp/C Y\ [-aj,bj], supple Y[ K\ (-bj,a,j) (3.8) 
ieN d jeN d 

i/ien i/ie Bedrosian identity \2. 1 0\) holds true. 

Proof: It suffices to prove that the integrand in each integral of (|3,5p vanishes almost everywhere. Let 
£ G Qk and j ^ k G N 2 d. We need to show that /(£ — i])g(rj) is zero for almost every rj £ Qj. Assume 
that there exists n G Qj that is in the support of /(£ — -)g. Then r\ G supp g n Qj and £ — 7/ G supp /. 
Since j 7^ fc, there exists I G such that ^rji < 0. We may assume that ^ > and r]i < 0. By (|3.8|1 . 
r// < —6/ while £1 — T}i < 6/, which is impossible. The contradiction completes the proof. □ 

To present the next sufficient condition, we set for j G N 2 d, := {A; G N 2 d : mj. = mj} and 
Qj := UkeijQk- We also call a subset AC1 11 closed under addition if for all x, y G A, x + y G A. 

Proposition 3.4 Let j G N 2 d. If supp / U supp (7 C Qj and Qj is closed under addition then 
f,g € L 2 ^^) satisfy the Bedrosian identity Ii2.10\) . 

Proof: By Theorem 13. 8} it suffices to show that for all £ G R d 

[ m-v)9(ri)(m(0-mj)dr ] = 0. (3.9) 
is, 

Clearly, the above equation holds true for £ G Qj. Assume that there exists £ G R d \ Qj that does not 
satisfy (13. 9p . Thus, there must exist 77 G Qj such that £ — 77 G supp / C Qj. Since Qj is closed under 
addition, £ — 77 G Qj and 77 G Qj imply that £ G Qj, a contradiction. □ 

In the one-dimensional case, the Bedrosian theorem has an appealing physical interpretation. 
Namely, it states that if / G L 2 (R) has low Fourier frequency and g G L 2 (R) has high Fourier 
frequency then they satisfy the Bedrosian identity fjl ,5j) . Similarly, condition (j3.8j) can be interpreted 
as that / has low Fourier frequency in each coordinate while g is of high Fourier frequency in each 
coordinate. It had been conjectured that the condition in the Bedrosian theorem was necessary for 
the one-dimensional Bedrosian identity (jl.5p until an explicit contradicting example was constructed 
in |36j. We recall that 

1 1 — 2t 2 
f(t) := -, g(t) := t G R 

satisfy H(fg) = fHg while supp / = supp g = R. Using this example, we define 

F(x) := J] f( Xj ), G{x) := ]J g( Xj ), x G R d . 

j£N d jgN d 

Since T is of the form (|2.2p . T(FG) = F(TG). This together with supp-F = suppG* = R rf implies 
that the condition (I3.8P is also unnecessary for the multidimensional Bedrosian identity (I2.10p . 

Surprisingly, a necessity about the one-dimensional Bedrosian theorem was obtained in [35] ■ It 
asserts that if /, g G L 2 (R) satisfy the Bedrosian identity (11.5p . supp / C [a, b] for some a,b > 0, and 
endpoints a, 6 are indeed contained in supp / then supp g CR \ (—6, a). We shall present an extension 
of this result for the multidimensional identity (|2.10p . To this end, we recall that the convolution ip*if; 
of ip,ip G L 2 (R d ) is defined by 

(tp* i/))(x) := / <p(x - t)ip(t)dt, x G R d . 
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It is well-known that suppc/? C supp ip + supp^. In the case that both tp,ijj are compactly 
supported, we have the celebrated Titchmarsh convolution theorem |17l 127] . 
For each subset A C R d , we denote by conv ^4 the convex hull of A in M. d . 



Lemma 3.5 If both ip,ip G L 2 (M. d ) are compactly supported then 

conv supp (if * ip) = conv suppy? + conv supp-0- 

Theorem 3.6 Let v be an extreme point of [-1, l] d , Q := {£ G M. d : > 0, j G N d }, Q' := -Q, 
m|g / wi|<3' ; and f,g G L 2 (]R rf ) satisfy supp/, suppg C Q U Q'. // a := (aj > : j G N^) and 
b := (6j > : j G N d ) that (ajvj : j G N d ), (—6^^ : j G N d ) G supp/ and 

supp / C {£ G Q : < a i; j G N d } U G -Q : < fy, j G N d }, 

i/ien ^ G L 2 (M. d ) satisfies the Bedrosian identity \2. 1 0\) if and only if 

supple {e^E^^ui^-^E^^j- 

Before moving on to the proof, let us understand the conditions in Theorem 13.61 Suppose that 
v = v 1 . Consequently, Q = Q\ and Q' = Theorem 13.61 states that given f,g G L 2 (IR rf ) with 
supp g C Qi U Q 2 , if 

supp / C [0, 0i ] |J Yl H»i,0], for some a = (a,,- > : j G N d ), 6 = (&,• > : j G N d ) 

jen d jen d 

and a, —b are actually contained in supp / then g satisfies the Bedrosian identity if and only if any 
£ G supp g lies in 

We illustrate the supports of / and g for the two-dimensional case in the following graph. 



t 2 



■ bt 



* t * * * * 



(a) supp/ (b) supp g 

We next present the proof of Theorem 13.61 
Proof: Without loss of generality, assume that Q = Qi and Q' = Q 2 - By the assumptions that 
supp/, supp g C Qi U Q2 and mi 7^ ni2, we obtain from Theorem 13.21 that g satisfies the Bedrosian 
identity (|2.10j) if and only if 

/ f(A-vMv)dv = o, C^Qi (3.10) 

JQi 
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and 



I f(Z-v)g(v)dv = o, £eQ 2 . 

JQa 



We shall only prove that supp (gXQ 2 ) — {£ e Qz '■ SjeN d 7T~ — — 1} smce the other inclusion that 

supp(<?x<2i) ^ {£, G Qi : Y2jm d I 1 — 1} can ^ e handled in a similar way. Set := fXQi an d 
-0 := gXQi- We decompose ^ uniquely into 

ip = ipi + i>2 

where suppV'i Q Y\j e n d [~ a j^] an d suppV>2 Q Q2 \ n,-eN d ( _a j> 0). By equation (|3.10p , 

(<p * il>i)(£) = -(<£* ^)(0> £ e Qi- 

Since 

supp (99 * C supp if + supp ^ 2 C Q 2 \ Yl [-aj, 0] + J| [0, a,] C| d \ Q°, 

jeN d jeN d 

where Q° is the interior of Q±, we have that 

suppO*Vl) CM d \Q°. 

This together with 

supp (tp C supp 99 + supp ?/> C J^J [— aj, aj] 

jeN d 

implies that 

(V9 * ^l) C {£ : ~aj < & < a,-, j G N d , ]T ^ < d - 1}. (3.11) 



conv 



jeN d 



We claim that for any £ £ supp^i, X^jeN d — ~~ 1- Assume to the contrary that this is invalid. Then 

CO 



there exists £° G conv supp^i such that X^jgN d IT > ~~ ^ Note that := a G supp(£>, which satisfies 
^ieN d 5~ = d - B y Lemma ESI 

conv supp (ip*ip{) = conv supp^i + conv supp(/9. 
It follows that there exists some point £ G conv supp (ip * tpi) such that 

6 



ieN d 3 



contradicting (|3.1ip . We conclude that 



suppV'i = supp (gXQ 2 ) ^ {£ G Q 2 : ~ - 

ieN d i 



which completes the proof. □ 

By a direct application of the Titchmarsh convolution theorem, one obtains another necessity 
theorem analog to the one-dimensional one proved in |35| . 
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Theorem 3.7 Let f G L 2 (R d ) satisfy for some a := (a,- > : j G N^) and 6 := (bj > : j G N^) i/iai 



supp / C [-aj,6j 



«7) 



d aZZ 7j/ie extreme points of Y[jeN d i~ a i> fy] are contained in supp/. T/ien (7 G £ 2 (M rf ) wi/j/i g being 



compactly supported satisfies the Bedrosian identity T(fg) = fTg for all operators T of the form \2. 2\) 
if and only if 

supp^C H R\(-bj, aj ). (3.12) 

jen d 

Proof: The sufficiency has been proved in Proposition 13.31 One sees that g satisfies T(fg) = fTg for 
all operators T given by (|2.2p if and only if for all j G 

Hj(fg) = fHj(g). (3.13) 



Set R d + := {£ G M d : > 0} and := {£ G R d : ^ < 0}. Introduce for each ip G L 2 {R d ) the 
associated pair of functions defined by 

3 \ 0, elsewhere, 5 \ 0, elsewhere, 
By Theorem 13.21 identity (|3.13p implies 

supp (f j+ * gjJ) C supp (/_,■_ * C R d j+ . 

By Lemma |3.5| we get 

supp (fj+) + supp (gj-) C supp (/,•_) + supp C 

Since the exterme points of rijGN d [ — a j-> bj] ne m supp/, there exist £ G supp/ J+ and n G /j_ with 
£j = 6j and r\j = —aj. This together with the above equation imply 

supp g C G M d : G R \ (-fy, a,). 

Since this is true for each j G N^, we get (j3. 12|) . □ 



4 Basic Multidimensional Analytic Signals 

Assume that an operator T of the form (|2.2p has been chosen to define multidimensional analytic 
signals. Applying T directly to a given signal would generally not yield physically meaningful results 
for the reason that the signal may contain multiple components. It is desirable to first decompose 
the signal into a sum of basic signals that behave well under T, and then apply the operator T to 
each summand. Note that T given by (|2.2p is a linear combination of compositions of partial Hilbert 
transforms. Thus, the purpose of this section is to construct basic signals that will behave well under 
each partial Hilbert transforms. More precisely, we shall first construct one-dimensional signals in 

M := {pcos9 : p G L 2 (R), 9 G C^K), p > 0, 9' > 0, H{pcos9) = psm9}. 

Multidimensional basic signals can then be formed by tensor products of functions in A4. 
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Constructions of functions pcosO in Ai has recently been considered in [23], where the phase 
function 8 was selected to be a strictly increasing function on M. Consequently, the constructed signal 
pcosO does not possess much fluctuation, which, on the other hand, is usually required in the time- 
frequency analysis. Thus, we shall adopt a different approach by choosing periodic phase functions 
that still enjoy nonnegative derivative. Specifically, we shall use phase functions determined by a finite 
Blaschke product. 

Set U := {z G C : \z\ < 1} and let a := (aj : j G N„) G U n . The finite Blaschke product B a 
associated with a is a holomorphic function on U defined as 

B a (z) := J] z G U. 



jeN n 

We let 6 a be the phase function determined by B a as 

e tfa(t) = B a ( e *), t G R, 

and aim at characterizing all real functions / G L 2 (M) such that 

H(f cos a ) = /sin0 a . (4.1) 

We start with a straightforward observation. 

Lemma 4.1 A real function f G L 2 (IR) satisfies 1\ ) if and only if 

H(fe i0 *) = -ife ie *. (4.2) 

Proof: Let / G L 2 (M) be real. Assume that (|4.ip holds true then we apply H to both sides of (|4.ip to 
get that 

-/cos (9 a = ff(/sin0 a ). 
Combining the above equation with (|4.1[) yields that 

F(/ e iea ) = ff(/cos0 a ) + zF(/sin0 a ) = /sin0 a -i/cos0 a = -i/e^, 

which is (|4.2p . Conversely, if (|4.2|) holds true then by comparing the real part of its both sides, we 
obtain (|4. 1|) . Thus, it suffices to show that (|4.2p is equivalent to (|4.3p . □ 

We first deal with the situation when a is a singleton {a}, a G U. In this case, we abbreviate 9 a 
as # a . Denote by r the backshift operator defined by 

(rf)(t) = f(t-l), teR. 

Lemma 4.2 If a consists of a single point a for some a G U i/ien a function f G L 2 (R) satisfies 
if and only if 

(rf)(0 = af(0, for all £ < 0. (4.3) 
Proo/: Let / G L 2 (R). We compute that 

(/e*»r = "«/ + (1 - l«| 2 ) E « fe_1 ^/. (4-4) 

fceN 

By taking the Fourier transform of both sides of (|4.2j) and engaging (jl.3j) . we see that (14. 2 p is equivalent 
to 
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By (|4.4p . the above equation can be rewritten as 

- o/(0 + (1 - |a| 2 ) ^ a fc ~ V/)(0 = 0, ? < 0. (4.5) 

fceN 

It remains to prove that (|4.5f) and (|4.3I) imply each other. Assume that (14.31) holds true. Then we 
observe for all £ < that 



-a/(0 + (1 - |a| 2 ) ^ a^V/tO = -a/(0 + (1 - |a| 2 ) ^ a*" 1 ^/® 

keN fcgN 

«/(0 ( -1 + (1 - |a| 2 ) XT 
V ken 



alY" 1 =0. 



On the other hand, assume that (|4.5|) holds. We apply r to the left hand side of (|4.5j) and then 
multiply it by a to obtain that 



3& 



|a| 2 (r/)(0 + (1 - |a| 2 ) £ a^V/XO = 0, £ < 0. 



fc=2 



Subtracting the above equation from (|4.5|) gives 

-a/(0 + |a| 2 (r/)(e) + (1 - |a| 2 )(r/)(£) = 0, £ < 0. 
An rearrangement of the above equation yields (14. 3p and completes the proof. □ 



We next present the crucial lemma leading to the characterization of real / satisfying (|4.ip . To 
this end, we denote for a := (aj : j £ N n ) 6 U ra by T a the operator 



(r a /) = ( II ( T - °S)) /■ 



v jeN„ 

Lemma 4.3 A function f G L 2 (M) satisfies {4.2fy if and only if 

(r a /)(e) = 0, /oraZZ£<0. (4.6) 

Proof: We shall use the induction on the number n of factors in the Blaschke product B a . By Lemma 
14.21 the result holds true when n = 1. Assume that n > 2 and the result also holds for the n — 1 case. 
Let a' := (a,- : 1 < j < n - 1) € U n_1 and 5 = fe ida ". Then (B~2j) is equivalent to 

H(ge i6 *') = -ige i9 *'. 

By induction, the above equation holds true if and only if 

(Ta'£)(f) = 0, for all £ < 0. 

We have computed that 

9 = -aj+(l - Kl 2 ) J^a^-V*/- 

fcGN 

Therefore, it suffices to show that fj4.6j) is equivalent to 

-a n T a >f(0 + (1 - K| 2 )tv ^a^ fc - 1 r fe /(£) = 0, for all £ < 0. 

fceN 
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Similar arguments as those in the proof of Lemma 14.21 fulfils this purpose. □ 

We conclude the results of this section by far into the following theorem. Note that i := r _1 is the 
shift operator defined as 

(Lf)(t) = f(t+1), t€R. 
We introduce another operator by setting 

jeN n 

Theorem 4.4 A real function f G L 2 (M) satisfies if and only if 

(rj)(0 = 0, for aU£<0 and (t-J)(0 = 0, for all £ > 0. (4.7) 
Proof: The result follow immediately from Lemmas 14.11 14.21 14.31 and the fact that for a real function 

/eL 2 (R),7(0 = /(-£UeK. D 

We next seek a closed form for real functions / satisfying (|4.ip . 

Theorem 4.5 Let Xj be the distinct elements in a with multiplicity rij, 1 < j < k. Then a real 
function f G L 2 (R) satisfies |^.i[ ) z/ and onZ?/ «/ i/iere exist functions Cji G L 2 (M) with supp(cjj)* C 
[—1,0], 1 < j < k, 1 < I < rij such that 



Proof: Assume that / is a real function in L 2 (R) that satisfies (|4.ip . Then by Theorem 14.41 / satisfies 
(|4.7|) . It implies that for each £ G (—1,0], /(£ — m) satisfies the difference equation 

(r a /)(C-m) = 0, mGZ + . (4.9) 

By the general solution of a difference equation, there exist functions eji on (—1,0], j G Z G N nj . 
such that 

l-X 



jeN k ZGN„ . q=l 



/(£ ~m)=J2 E e i* W [I " +^)' m G Z +- ( 4 - 10 ) 
Since / is real, we have that 

f(x) =-L= [ f^) e ^ = -Re( f f(0e**d£ 

V27T Jm 7T \J-oo J 



^ Re (E / KOe^dA (4.11) 

-Re E / /(e-m^^de 
\n=cr 1 



Now we apply ([4.10|) to obtain that 

i-i 



E/ /(C-m) e -«-)^=E E / E A " e_ima; II( m + 5)) e iK0^. (4.12) 

mMK -1 jeN fc ZgN„ V -1 ^m=0 g=l ' 



15 



Recall that 

oo 1-1 (7 -IV 

£ x ? e ~ mx IT " + ?) = a : A . e -^v J e Nfc ' 1 e Nn *- (413) 

m=0 9=1 \ 3 J 

Combining (liTTT) . (12T21) . and (|4TT51) yields that 
Finally, setting 

Cil (x) = / e^COe^de, i E N fcl Z E N n , 

leads to the form (|4.7p of /, 

Conversely, if / os of the form (|4.7p then by reversing the above arguments, we obtain (|4,10p . 
Thus, for each £ £ (—1, 0], /(£ — m) satisfies the difference equation (|4.9p . Therefore, / satisfies (|4.6p . 
Since / is real, (|4.7p holds true. By Lemma 14.41 / satisfies (|4.1|) . □ 
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